 use "Data/DawesBond_European_AAO2023(RFS).dta", clear

 keep date P_*_GBP
 
 
gen int year = yofd(date)
 
 
collapse (mean)  P_london_GBP_GBP P_paris_GBP_GBP P_amsterdam_GBP_GBP P_zurich_GBP_GBP, by(year)

save "Temp/bonds_annual", replace 


 import excel "Data/DawesBond_PS2023(JPE).xlsx", sheet("Sheet1") firstrow clear

 tsset Datesweekendson
 
 gen int year = yofd(Datesweekendson)

 collapse (mean) Dawes_NYSE, by(year)
 
 
 merge 1:1 year using "Temp/bonds_annual"
 
 
 keep if year==1936
 
 rename Dawes_NYSE P_NYC
 rename P_london_GBP_GBP 		P_London
 rename P_paris_GBP_GBP  		P_Paris
 rename P_amsterdam_GBP_GBP    P_Amstersdam
 rename P_zurich_GBP_GBP		P_Zurich
 
 
 reshape long P_, i(year) j(market) string
 
 
 
 
 
 gen barorder=1 if market=="London"
 replace barorder=2  if  market=="NYC" 
  replace barorder=3  if  market=="Paris" 
 replace barorder=4  if  market=="Amsterdam" 
 replace barorder=5  if  market=="Zurich" 
 
 replace market="New York" if market=="NYC"
 
 
 
 graph bar P_  , over(market, sort (barorder) ) ytitle("Bond price in % of par (1936)") label xsize(3) ysize(4)
 
 
 graph export "Output/Figure2a.pdf", replace 

 
 
 
 